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ABSTRACT 

X-ray spectroscopic measurements of the Cygnus Loop supernova remnant 
indicate that metal abundances throughout most of the remnant's rim are de- 
pleted to ~0.2 times the solar value. However, recent X-ray studies have revealed 
in some narrow regions along the outermost rim anomalously "enhanced" abun- 
dances (up to ~1 solar). The reason for these anomalous abundances is not un- 
derstood. Here, we examine X-ray spectra in annular sectors covering nearly the 
entire rim of the Cygnus Loop using Suzaku (21 pointings) and XMM-Newton (1 
pointing). We find that spectra in the "enhanced" abundance regions commonly 
show a strong emission feature at ~0.7keV. This feature is likely a complex of 
He-like O K(7 + 6 + e), although other possibilities cannot be fully excluded. 
The intensity of this emission relative to He-like O Ka appears to be too high to 
be explained as thermal emission. This fact, as well as the spatial concentration 
of the anomalous abundances in the outermost rim, leads us to propose an origin 
from charge-exchange processes between neutrals and H-like O. We show that the 
presence of charge-exchange emission could lead to the inference of apparently 
"enhanced" metal abundances using pure thermal emission models. Account- 
ing for charge-exchange emission, the actual abundances could be uniformly low 
throughout the rim. The overall abundance depletion remains an open question. 
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1. Introduction 

The importance of charge- exchange (CX) emission in X-ray astrophysics emerged with 
the discovery of bright X-ray emission from comet Hyakutake (Lisse et ah 1996) and the sub- 
sequent identification of its origin from CX processes between neutral atoms in the comet's 
atmosphere and highly charged ions in the solar wind (Cravens et al. 1997). A number of 
comets are now known to have X-ray emission originating from CX processes (e.g., Cravens 
2002). More recently, the diffuse soft X-ray background, which was previously thought to 
be dominated by thermal emission from the local hot bubble, turned out to be significantly 
contaminated by CX emission from interaction of heliospheric/geocoronal neutrals and solar 
wind ions (e.g.. Cox 1998; Cravens 2000; Wargelin et al. 2004; Snowden et al. 2004; Lalle- 
ment 2004a; Fujimoto et al. 2007; Smith et al. 2007; Koutroumpa et al. 2009; Ezoe et al. 
2010). 

In principle, CX-induced X-ray emission could be produced at any astrophysical site 
where hot plasma interacts with (partially) neutral gas. One very promising site is the thin 
post-shock layer in supernova remnants (SNRs), since SNR shocks are collisionless and both 
unshocked cold neutrals and shocked hot ions can be present just behind the shock front. 
In fact, observational evidence of CX emission has been obtained as the broad component 
of Ha emission in many SNRs for more than 30 yrs (e.g.. Chevalier et al. 1980; Ghavamian 
et al. 2001). On the other hand, CX-induced "X-ray" emission has not yet been detected 
firmly. It has been suggested only for the SMC SNR 1E0102.2-7219 (E0102) to explain 
anomalously high (in the context of thermal emission) intensity ratios of and 7//? in 
O VIII Ly series as observed by the XMM-Newton Reflection Grating Spectrometer (RGS) 
(Rasmussen et al. 2001). The authors suggested that CX processes at the reverse shock 
contribute significantly to the X-ray emission. However, Chandra High Energy Transmission 
Grating (HETG) observations (Flanagan et al. 2004) could not confirm the RGS result: 
the /3/a ratio was found to be about two thirds of that measured by the RGS, consistent 
with that expected in thermal emission with an electron temperature of ~0.5 keV. Although 
the detection of CX-induced X-ray emission from E0102 has been debated, other SNRs, 
especially nearby Galactic remnants, are more suitable for detecting it. Taking account of 
the CX emission could be extremely important for understanding the true properties of the 
shocked plasma, since the undetected presence of CX emission leads to incorrect plasma 
parameters (e.g., the metal abundance and the electron temperature) if we interpret the 
CX-contaminated X-ray spectra with pure thermal emission models. It should be also noted 
that CX reactions may be important for understanding magnetic-field amplification behind 
the shock front, a topic of continuing interest in SNR shock physics (Ohira et al. 2009; Ohira 
& Takahara 2010). 
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Unusual metal-abundance regions have been reported at the rim of the Cygnus Loop, 
which is a nearby (540 pc: Blair et al. 2005), extended (diameter of ~ 2°. 8: Levenson et 
al. 1998) SNR: the metal abundances in some of the outermost narrow (a few arcminutes, 
or <0.5pc at a distance of 540 pc) rim regions are anomalously "enhanced" (up to about 1 
solar, hereafter we define these abundances as "enhanced", Katsuda et al. 2008ab; Tsunemi 
et al. 2009; Uchida et al. 2009; Kosugi et al. 2010), whereas the metal abundances are 
typically around 0.2 times the solar value for most of the rim (hereafter we define these 
abundances as "normal" , Miyata et al. 1999; Leahy 2004; Tsunemi et al. 2007; Miyata et al. 
2007; Nemes et al. 2008; Zhou et al. 2010). Above and throughout this paper, we use the 
solar abundances of Anders & Grevesse (1989). The abundance ratios among metals in the 
"enhanced" abundance regions are consistent with solar values within a factor of ~2. This 
fact, as well as the localization of the "enhancements" to a few regions along outermost rim, 
led the authors to suggest that the "enhancements" arc not likely due to contamination by 
SN cjecta (Katsuda et al. 2008a). It is believed that the Cygnus Loop is a remnant from a 
corc-coUapse SN and that its forward shock is now hitting the wall of a wind-blown cavity over 
considerable fraction of the rim (e.g., McCray & Snow 1979; Hester et al. 1994; Levenson et 
al. 1997). Therefore, stellar winds might have altered abundances in part of the surroundings. 
However, stellar winds are expected to be rich in only fighter elements such as C, N, O, and 
Ne, as a result of CNO processing and/or He burning (e.g., Rauscher et al. 2002; Murashima 
et al. 2006), which is inconsistent with the overall abundance "enhancements". It is also 
unlikely that the ISM abundances are inhomogeneous over such a small scale (a length scale 
of several pc), given that ISM abundances have been reported to be fairly constant in the 
local Galaxy. For example, Cartledge et al. (2004; 2006) revealed that the ISM abundances 
of a number of sight lines shorter than 800 pc are well within the measurement uncertainties, 
and concluded that the intrinsic rms scatter of the ISM abundances is ~10%. A similar 
rms scatter of the ISM abundances was recently obtained from observations of nearby B 
stars (Przybilla et al. 2008). Since the Cygnus Loop, whose distance is 540l8Q°pc (Blair et 
al. 2005), is within the region investigated, the abundance variation around the rim of the 
Cygnus Loop is difficult to understand. In addition to this problem, the "normal" abundances 
(~0.2 times the solar value) are also puzzling, given that the ISM metal abundances around 
the Cygnus Loop are about half the solar value (e.g., Cartledge et al. 2004). 

Noting that the "enhanced" abundances are close to the ISM abundances, Katsuda et al. 
(2008b) and Tsunemi et al. (2009) suggested that the "normal" (i.e., depleted) abundances 
are only apparently low due to presence of a nonthermal component (power-law continuum), 
which can artificially raise the continuum level and thus decrease the equivalent widths of the 
lines. Indeed, introduction of a nonthermal component does increase the fitted abundances 
and yields sfightly better overall spectral fits (the F-test probability of 99%). However, 
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the radio flux inferred by tlie X-ray emission is unrealistic: it is required to be at least a 
few orders of magnitude higher than that observed. As a result, the authors speculated 
that there are indeed two kinds of material with different abundances: one is the ISM 
with "enhanced" abundances and the other, having the "normal" abundances, is the wall 
of the cavity surrounding the SN precursor. There is no good reason, however, why the 
ISM abundances should be elevated compared with those of the cavity wall, casting some 
doubt on this interpretation. Although some other possible explanations of the abundance 
enhancements and/or depletion, including dust depletion and resonance-line scattering, have 
been discussed in the literature (e.g., Katsuda et al. 2008b; Miyata et al. 2008; Uchida et al. 
2009; Kosugi et al. 2010), none provides a satisfactory explanation. 

In this paper, we suggest that the presence of CX emission is responsible for the abun- 
dance "enhancements" in the outermost rims of the Cygnus Loop. We present a new 
spatially-resolved X-ray spectral analysis of nearly the entire rim of the Cygnus Loop using 
Suzaku (21 pointings) and XMM-Newton (1 pointing). We argue that CX- induced X-ray 
emission is significant in the "enhanced" abundance regions and not in the "normal" abun- 
dance regions, and that it introduces apparent abundance enhancements when we apply pure 
thermal emission models to the CX-contaminated spectra. 

2. Observations and Data Reduction 

We have observed nearly the entire rim region of the Cygnus Loop with Suzaku. In 
our analysis, we concentrate on the XIS data, since hard X-ray emission (>10keV) from 
the Cygnus Loop is negligible. The fields of view (FOV) of the XIS pointings are shown as 
white boxes overlaid on the X-ray mosaic of the Cygnus Loop obtained by the ROSAT HRI 
(Fig. [1]). We also analyze the XMM-Newton EPIC data. XMM-Newton' s FOV, which is 
indicated as a white circle in Fig. [T|, partly covers the gap in the Suzaku observations along 
the eastern rim. Information about the observations is summarized in Table [H 

We reprocessed the XIS data using the 20090615 version of the CTI calibration file. We 
filtered the reprocessed data using the standard criterig0 recommended by the Suzaku/Xl'^ 
calibration team. As for XMM-Newton, the raw data were processed using version 8.0.0 
of the XMM Science Analysis Software (SAS). We selected X-ray events corresponding to 
patterns 0-12@, and removed all the events in bad columns listed in Kirsch (2006). After the 



"'^See the Suzaku data reduction manual which can be found at 
http: / /heasarc . gsfc . nasa. gov /docs /suzaku / analysis / abc[ 

^http://heasarc.nasa.gov/docs/xmni/abc/ 
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filtering, the data were vignetting-corrected using the SAS task evigweight. The effective 
exposure times for the various observations are hsted in Table [1] 

Background in the XIS arises from two sources: non- X-ray background (NXB) caused 
by charged particles and 7-rays hitting the detectors, and diffuse X-ray background from the 
sky. We generate NXB spectra suitable for our observations using the xisnxbgen software 
(Tawa et al. 2008), and subtract them from the source spectra. The X-ray background was 
extensively studied by Yoshino et al. (2009) who modeled it as two thin thermal components 
plus one (or two) power-law component (s). They examined 15 blank-sky regions in various 
directions. Among these, the closest one to the Cygnus Loop located at (^, 6) = (74°.0, -8°. 5) 
is "Low Latitude 86-21" located at {i, 6) = (86°.0, -20°. 8). In our analysis, we model the 
X-ray background based on the spectral parameters for this blank-sky region. We note that 
the X-ray background is much weaker than the Loop emission (which will be shown in e.g.. 
Fig. |3]) and that the ROSAT all-sky survey shows generally similar backgrounds in these 
regions (Snowden, et al. 1997). Therefore, any differences in the background are negligibly 
small. The parameters of the background are taken from Table 4 in Yoshino et al. (2009) 
except for the total Galactic hydrogen column density which we modified to 1.7xl0^^cm~^ for 
the direction of the Cygnus Loop (Leiden/Argentine/Bonn Survey of Galactic HI, Kalberla 
et al. 2005). For the EPIC background, we use the data set accumulated from blank sky 
observations prepared by Read & Ponman (2003). We subtract the blank-sky data from the 
source after matching the detector coordinates. 

3. Comparison of "Enhanced" and "Normal" Abundance Regions 

As described in the introduction, there are apparent abundance inhomogeneities along 
the rim of the Cygnus Loop: overall, the abundance is typically 0.2 times the solar value, 
while parts of the outermost rim show "enhanced" abundances of up to ~1 solar. Our inves- 
tigation of the nature of the enhancement starts by summarizing observational properties of 
the rim regions. 

3.1. Spectral Properties 

First, we compare a typical "enhanced" abundance spectrum with a "normal" one. 
To this end, we focus on the XIS data for the NE rim, because these were obtained in 
the early phase of the Suzaku mission, and have the best photon statistics as well as the 
best spectral resolution among all the data listed in Table [TJ These data were previously 
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analyzed by Katsuda et al. (2008a), who divided each FOV into small box regions, and found 
"enhanced" abundances in the outermost regions in NE3 and NE4. They also found that the 
parameters are generally constant with azimuth (along the shock front). Therefore, to make 
the comparison with better statistics, we here extract azimuthally-integrated spectra from 
narrow regions (thickness of 2') along the shock front in each FOV as shown in Fig.[Tl Figure[2] 
shows spectra extracted from the outermost regions in NE4 (for the "enhanced" abundance 
region) and NE2 (for the "normal" abundance region). These spectra are normalized to 
each other in the energy bands of 1.4-1.5 keV and 1.65-1.75 keV where line features are not 
evident. The differences between the two are readily seen: K-shell lines from O and Ne are 
significantly enhanced in NE4 relative to those in NE2, whereas the complex of Fe L-shell 
lines at ~0.82keV is comparable in the two spectra. 

To confirm previous fitting results, we fit the two spectra with an absorbed, single 
component, plane-parallel shock model [a combination of the Tbabs (Wilms et al. 2000) and 
the vpshock model (NEI version 2.0) (e.g., Borkowski et al. 2001) in XSPEC v 12.5.1]. This 
thermal emission model is the same as that employed in Katsuda et al. (2008b). We assume 
the hydrogen column density of the intervening material to be 3x10^° cm^^ (e.g., Kosugi et 
al. 2010). Note that this value is much smaller than the total Galactic column density of 
1.7xl0^^cm~^ for the direction of the Cygnus Loop (see. Section 2) because of the proximity 
of the Loop. We freely vary the electron temperature, kTe, the ionization timescale, riet, and 
the volume emission measure (VEM = J nenndV, where rie and tt-h are number densities 
of electrons and protons, respectively, and V is the X-ray-emitting volume). Above, Ugt 
is the electron density times the elapsed time after shock heating, and the model uses a 
range from up to a fitted maximum value. The abundances of several elements, whose line 
emission is prominent in the spectrum, are treated as free parameters: C, N, O, Ne, Mg, Si 
(=S), Fe (=Ni). We manually adjust the energy scale within the systematic uncertainties 
to obtain better fits, by allowing energy-scale offsets to vary freely for the front illuminated 
(FI) and the back illuminated (BI) detectors, respectively. In this way, we fit the spectra 
in energy ranges of 0.4-3.0 keV for FI and 0.33-3.0 keV for BI. Figure |3] shows the best- 
fit model along with the data points. The best-fit parameters are summarized in Table |2j 
We confirm that the absolute (relative to H) abundances are enhanced in NE4 compared 
with those in NE2, and that the abundance ratios of C/Fe, N/Fe, 0/Fe, and Ne/Fe are 
~2-4 times higher in NE4 than those in NE2. We also confirm that the plasma is in the 
non-equilibrium ionization (NEI) condition {net < 10^^ cm~^ sec~^). Dividing the ionization 
timescale by the time elapsed after shock heating, which is roughly inferred from the distance 
of ~ 1 X 10^^ cm (i.e., 2' at a distance of 540 pc) divided by the fluid velocity of ~200 kmsec"^ 
(e.g., Salvesen et al. 2009), we estimate the post-shock electron densities in NE4 and NE2 
to be ~1 cm~^ and ~6cm~^, respectively. Assuming strong shocks, we can estimate the 
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pre-shock ambient densities in NE4 and NE2 to be ~0.25cm~^ and ~1.5cm~^, respectively. 
These values generally agree with previous estimates (e.g., Raymond et al. 2003 for NE4; 
Hester et al. 1994 for NE2), showing the robustness of the ionization timescales derived in 
these outermost rim regions from the very shock fronts to 2' inside the shocks. 

We notice that the fit quality for the NE4 spectrum is far from formally acceptable 
(reduced of ~2.2), while the fit quality for the NE2 spectrum is much better. The NE4 
spectrum in Fig. [3] (a) shows apparent excess emission at ~0.7keV relative to the best- 
fit model. This discrepancy was already noticed in the Chandra analysis of the NE rim 
(Katsuda et al. 2008b), although it was not so evident due to poorer statistics. We therefore 
re-fit the NE4 spectrum using the same model, but ignoring the energy band 0.68-0.76 keV. 
Table [2] (the 3rd column) summarizes the best-fit parameters. The fit quality is significantly 
improved over the original one (the reduced-x^ is reduced from 2.2 to 1.6). Both the electron 
temperature and the abundances significantly decrease from the original best-fit values to 
values more consistent with those in NE2 (the 2nd column in Table [2]). These changes 
are interpreted as follows. In the original (entire energy band) fitting, the He-like O K- 
shell lines (n > 4 — )■ n = 1) which dominate around 0.7keV are forced to be as strong 
as possible to match excess emission in the 0.68-0.76 keV band. The way these lines are 
made stronger is by an increase in the electron temperature as well as the O abundance. 
When this energy band is disregarded in the fit, the electron temperature decreases from 
~0.3keV to ~0.2keV. The reduced temperature in turn allows a reduced O abundance, 
since the emissivities of O K-shell lines peak at kT^ ~0.2keV. Other metal abundances are 
also reduced, with the net result that the relative abundances among metals are maintained. 
Therefore, we conclude that the 0.68-0.76 keV band plays an important role in determining 
absolute abundances in the "enhanced" abundance spectrum when using the Suzaku XIS 
or any other CCD spectrometers. For the "normal" abundance spectrum (e.g., the NE2 
spectrum in Fig. [3] (b)), ignoring the 0.68-0.76 keV band does not change the fit quality or 
the best-fit parameters. 

In order to compare the two spectra in more detail, we next fit them with a phe- 
nomenological model consisting of an absorbed, bremsstrahlung continuum plus a number 
of Gaussians for line emission from C, N, O, Ne, Mg, Si, and Fe. Again, we assume the 
hydrogen column density to be 3x10^° cm~^. We include eight Gaussians for the prominent 
line features at -0.35 keV (C Lya), -0.45 keV (C Ly/3 + N Hea), -0.5 keV (N He^ + N 
Lya), -0.57keV (O Hea), -0.91 keV (Ne Hea), 1.35 keV (Mg Hea), -1.6 keV (Mg He/3), 
and — 1.85keV (Si Hea). The center energies of these lines are allowed to vary freely. In 
addition to these, we introduce nine Gaussians to represent lines that could contribute in 
the crowded 0.7-1.2 keV band. These are O Lya at 0.654keV, O He/3 at 0.666 keV, an Fe 
XVII L-shell complex from the 3s— 7'2p transitions at 0.726 keV, O Ly/3 at 0.775 keV, an Fe 
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XVII L-shell complex from the 3d^2p transitions at 0.822 keV, Ne Lya at 1.022 keV, Ne 
He/3 at 1.074keV, Ne He7 at 1.127keV, and Ne Ly/3 at 1.210keV. The hne center energies 
of these nine Gaussians are fixed at their theoretically expected values (APED: Smith et al. 
2001). Based on the most recent NEI plasma code (Borkowski et al. 2001; the augmented 
NEI version 2.0 that includes inner-shell processes), the intensity ratio between overall Fe 
L-shell 3s— )-2p and 3d— >-2p lines is fairly constant at unity. Therefore, we assume the normal- 
ization of the Gaussian for the 3s— )-2p transition lines to be the same as that for the 3d— )-2p 
transition lines. We manually adjust the energy scale within the systematic uncertainty. In 
this way, we fit the two spectra shown in Fig. [2l As expected from the spectral fitting in 
the previous paragraphs, we find excess emission from the best-fit model at ~0.7keV for the 
"enhanced" abundance spectrum. Addition of a Gaussian at ~0.7keV gives a satisfactory 
fit; the statistical significance of introducing the line is greater than 99.9% for both spectra 
(reduced-x^ values are reduced from 1297 to 601 for NE4 and from 440 to 404 for NE2). 

The electron temperatures of the bremsstrahlung components are 0.16±0.01keV and 
0.26±0.01keV for NE2 and NE4, respectively. Although we do not consider radiative re- 
combination continuum emission or two-photon decay continuum emission, the electron tem- 
peratures derived are reasonable for the rim of the Cygnus Loop. Table [3] summarizes details 
of this model. The line center energies listed in the table include the effects of energy-scale 
shifts determined by the fittings. Figure H] presents the best-fit models. We find that the 
line intensities in NE4 are generally ~3 times higher than those in NE2, whereas the line 
at ~0.7keV indicated as a red line in Fig. HJis ~7 times stronger in NE4 than in NE2. Its 
nonzero width suggests that it is a complex of lines, although we cannot fully eliminate the 
possibility that the broadening might be due to calibration effects. 



3.2. Spatial Properties 

Before proceeding to an investigation of the origin of the ~0.7keV feature, we briefly 
summarize the spatial properties of the ~0.7keV emission. To reveal the distribution of 
the ~0.7keV feature, we fit all the spectra for all the regions in Fig. [1] using the same 
phenomenological model as we employed in Section 3, a bremsstrahlung plus 18 Gaussians. 
We show example best-fit values for the outermost rim regions in P27, P19, Rim2, and P24 
in Table [31 After deriving individual line intensities in each region, we calculate the line ratio 
between the ~0.7keV emission and the Fe L-shell complex at ~0.82keV. This ratio should 
be a good tracer of the relative strength of the ~0.7keV emission, given that the abundance 
(and intensity) of Fe (L-shell complex) is relatively constant in the rim regions (see, e.g.. 
Fig. |2]). Figure |5] (a) shows this line-ratio map, with the X-ray boundary of the Cygnus Loop 
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in white. The map clearly shows that the ~0.7keV emission is enhanced exclusively within 
a narrow (< 4') layer behind the shock. 

Also notable from Fig.|5](a) is the azimuthal variation. The ~0.7keV emission is evident 
only at position angles of 0°-40° (N-NE), 110°-160° (SE), and 270°-330° (northwest: NW) 
measured from north over east. These are the same locations as where the "enhanced" 
abundances have been reported so far (N-NE by Katsuda et al. 2008ab and Uchida et al. 
2009; SE by Tsunemi et al. 2009 and Kosugi et al. 2010). In addition, preliminary spectral 
analysis of the NW region shows "enhanced" abundances at the outermost regions compared 
with those in the inner regions (Takakura et al. in preparation). Therefore, it appears that 
the "enhanced" abundance regions commonly show the strong ~0.7keV emission. We further 
note that a map of the intensity ratio of Ne Hea to the Fe L-shell complex at ~0.82keV 
is quite similar to Fig. E] (a), and that the intensity-ratio enhancements seen in part of the 
region (which corresponds to the "enhanced" abundance regions) are statistically significant. 

Figure [5] (b) shows a 325 MHz radio image of the Cygnus Loop from the The Westerbork 
Northern Sky Survey (WENSS) (Rengelink et al. 1997) covering the same sky region of Fig.O 
(a). Looking at Fig. \5\ (a) and (b), we notice that the regions of strong ~0.7keV emission 
are generally anti-correlated with radio-bright regions. 

4. Discussion 

4.1. Identification of Excess 0.7 keV Emission in the "Enhanced" Abundance 

Regions 

In this section, we discuss the origin of the ~0.7keV emission feature. The center 
energy of the ~0.7keV feature suggests that it is most likely either O K-shell emission or Fe 
XVII/XVI L-shell emission. We first consider Fe XVII. In the X-ray band, there are two well- 
known strong Fe L-shell complexes at ~0.73keV and ~0.82keV, which mostly arise from the 
3s— )-2p and the 3d— >2p transitions in Fe XVII ions, respectively. While these two complexes 
are included in the phenomenological model above, the intensities of the two were set equal, 
based on the most up-to-date plasma code (Borkowski et al. 2001). This assumption leaves 
room for the 3s— i'2p complex to contribute to the ~0.7keV complex, given the moderate 
spectral resolution of the XIS (FWHM~60eV at IkeV). However, the energy of the 3s— ^2p 
complex is ~0.73keV, making it unlikely to be the dominant source of the ~0.7keV emission, 
the center energy of which is 0.705±0.002stat i O.OOSgygkeV (the systematic uncertainty is 
from Koyama et al. 2007). In addition, if we attribute the entire intensity of the ~0.7keV 
emission to the 3s— >2p complex, the intensity ratio of 3s— >^2p/3d— >-2p would be ~3.5. Such a 
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high ratio has never been observed in an astrophysical object. Although astrophysical sources 
including stellar coronae, SNRs, and elliptical galaxies do tend to show higher 3s— )-2p/3d^2p 
ratios than the predicted value (Doron & Behar 2002 and references therein), the observed 
ratio ranges from 1 to 2. Our own analyses at XMM-Newton RGS observation of the LMC 
SNRs DEM L71, N103B, N132D, N49, and N63A reveal that the 3s^2p/3d^2p intensity 
ratio in each is near or less than unity. Therefore, we rule out 3s^2p emission from Fe XVII 
as the origin of the ~0.7keV feature. 

Turning to Fe XVI L-shell lines, Graf et al. (2009) recently reported relative intensities, 
using University of California / Lawrence Livermore National Laboratory's EBIT-L In the 
energy range of 0.69-0.86 keV, the authors measured ten Fe XVI lines originating from inner- 
shell ionization processes. Among them, the strongest one is located at ~0.82keV, whereas 
four relatively weak lines are clustered around 0.7 keV. The total intensity of the lines around 
0.7keV is weaker than that around 0.82 keV. This measurement is inconsistent with our mea- 
sured line ratio: the ~0.7keV emission is ~2.5 times stronger than the ~0.82keV emission 
(which should include both the Fe XVI L-shell lines and the Fe XVII 3d^2p complex). It 
is therefore unlikely that the Fe XVI L-shell lines dominate the ~0.7keV emission, although 
we cannot fully exclude this possibility given that these lines will incidentally be stronger in 
an NEI plasma than in a coUisional ionization equilibrium plasma. 

Finally, we consider the possibility that O VII K-shell lines (n > 4 — )■ n = 1) produce 
the ~0.7keV emission. In this case, we again run into a line ratio problem for a thermal 
spectrum: the intensity ratio of the ~0.7keV emission to O Hea is measured to be ~0.06 
(see, the parameters in Table [3]), whereas the O He(7 + 5) to O Hea ratio is expected to 
be ~0.02 at kTe = 0.2 keV based on the NEI code (Borkowski et al. 2001). Therefore, it 
is also unlikely that thermally-emitted O VII K-shell lines are responsible for the ~0.7keV 
emission. 



4.2. Charge-Exchange Emission 

We suggest here that a pure NEI plasma model is not appropriate for reproducing the 
~0.7keV emission in the "enhanced" abundance region, and that an additional emission 
mechanism is required. One possible emission mechanism, suggested from the interpretation 
of the XMM-Newton RGS spectrum of the SMC remnant E0102 (Rasmussen et al. 2001), is 
CX-induced emission between neutrals and ions at the forward shock of the Cygnus Loop. 

CX emission would play an important role only in a thin region immediately behind 
the forward shock in a SNR, because cold neutrals cannot survive far behind the shock. 
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Taking account of line-of-sight effects, CX emission if present could only be detected at the 
periphery of a SNR. This is exactly what the line ratio map in Fig. |5] (a) shows. 

In general, the presence of neutral hydrogen is essential for CX emission to be sig- 
nificant. The presence of neutral hydrogen in the post-shock region can be tested by the 
presence/absence of non-radiative Ha filaments. Such filaments are found around nearly the 
entire perimeter of the Cygnus Loop (Levenson et al. 1998), clearly indicating the presence 
of neutral hydrogen in the post-shock regions. One Ha filament along the NE rim was exten- 
sively studied by Ghavamian et al. (2001), who modeled the broad and narrow components 
of the Ha and H/3 emission with their shock model to measure the degree of electron-ion 
temperature equilibration. In their model, the neutral fraction in the ambient medium is 
also a measurable parameter, although not a sensitive one. They examined three cases of 
neutral fractions, i.e., 0.1, 0.5, and 0.9, to model the broad-to-narrow line ratios, and found 
that the first case is slightly outside the range of the observational uncertainties whereas 
the latter two are consistent with the observation. Thus, the preferred range for the neutral 
fraction is larger than 0.1 in this region. 

With non-radiative Ha filaments around nearly the entire periphery, and a non-zero 
neutral fraction, it might be expected that the CX emission should be closely correlated 
with the non-radiative filaments. However, it does not seem to be the case. One possible ex- 
planation for distribution of CX emission is suggested by the radio map (Fig. |5]). The regions 
of the rim from which the hypothetical CX emission is not detected generally correspond to 
those that are radio bright. The radio emission is likely produced by accelerated cosmic-rays 
(electrons), since the spectral shape of the radio emission is well described by a power-law 
continuum (e.g., Uyaniker et al. 2004). Therefore, the anti-correlation can be qualitatively 
understood in a view that the cosmic-ray shock precursor more effectively pre-ionizes the 
medium, and thus reduces the CX emission below a detectable level. It is also interesting to 
note that the radio-bright (or CX-absent) regions correspond to radiative filaments (Leven- 
son et al. 1998). This correspondence suggests that the ambient densities there are higher 
and the shock velocities are lower than elsewhere along the rim. According to Lallement 
(2004b), slower shocks produce less CX emission, providing another possible explanation for 
the observed azimuthal variation of the CX emission. It may be also possible that thermal 
emission overwhelms CX emission in these radio/optical-bright regions. 

All of these observations support the possibility of CX emission in the "enhanced" abun- 
dance regions. Furthermore, a major benefit of this scenario is that the introduction of CX 
emission resolves the mysterious abundance "enhancements" in the outermost rim regions, 
by reducing the "enhanced" abundances to consistency with the "normal" abundances found 
everywhere else around the rim. We therefore suggest that CX emission is contributing sig- 
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nificantly to the X-ray spectrum in the "enhanced" abundance regions. 

We note that if this suggestion is to hold, CX emission other than the ~0.7keV complex 
is expected to be present. As pointed out in Section 3.2, the similarity of a map of the 
intensity ratio of Ne Hea to the Fe L-shell complex to Fig. [5] suggests that Ne Hea is 
also contaminated by CX-induced lines. In addition, under the condition of the electron 
temperature of ~0.2keV, C is likely fully ionized and thus H-like C emission due to CX 
processes would be also significant. CX strongly enhances high principal number Ly emission 
such as C Ly7 at ~0.46keV (Wargelin et al. 2008). This emission can contaminate the C 
LyP+N Hea blend around 440 eV (see. Table [3]). Therefore, the line intensity ratio of this 
blend to C Lya may be enhanced where we find indications of CX. However, such a trend 
can not be detected as shown in Table [31 A definitive identification of the possible CX lines 
would require higher spectral resolution observations which will be soon available by the Soft 
X-ray Spectrometer (SXS) with the energy resolution of 4-7 eV and the energy calibration 
of 1-2 eV (Mitsuda et al. 2010) onboard the A stro-H satellite (Takahashi et al. 2010). 



4.3. Supporting Evidence for CX Emission 

4-3.1. Spectral Modeling of CX Emission 

In order to further study the properties of the proposed CX emission, we model CX- 
contaminated spectra by a combination of a number of Gaussians for the CX emission and 
a plane parallel shock model (the vpshock model; Borkowski et al. 2001) for the thermal 
emission. We examine three spectra extracted from the outermost rim regions in NE4, P27, 
and Rim2 where the possible CX emission is fairly strong (see. Fig. |5]). The challenge is 
the difficulty in separating CX emission from thermal emission in CCD spectra. Detailed 
modeling considering the effects of CX on the ionization balance or evolution of the electron 
temperature is far beyond the scope of this paper. We here crudely assume that the shape 
of the thermal emission in this region is similar to that in its surroundings, so that we can 
estimate the thermal emission in the CX-contaminated spectra. Since the thermal spectra 
are expected to be more constant with azimuth than with the distance from the shock front, 
we use the thermal spectrum extracted from the outermost rim in observation fields where 
we do not find indications of CX emission. As a template of thermal emission, we take the 
NE2 spectrum shown in Figs. [2] and [3l The best-fit spectral parameters are listed in Table [2l 
It should be noted that recombination rates of ions increase due to CX processes and thus 
plasmas showing CX emission may be in lower ionization states than those showing pure 
thermal emission. However, it is quite difficult to estimate this effect quantitatively. We 
thus evaluate systematic uncertainties by examining two ionization timescales: one is the 
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value derived in the NE2 spectrum and the other is half the value. 

In the fitting, we fix all the thermal parameters except for the emission measures (EM 
= / neUiidi, where i is the plasma depth). We introduce 22 Gaussians to represent possible 
CX line emission. These are K-shell lines from C, N, O, Ne, Mg, and Si. For simplicity 
of the model, we ignore Fe L-shell lines that might be significant in the X-ray band. For 
example, Fe L-shell lines may affect intensities of the O Ly(/3 + 7 + 5 + e) around 0.8 keV, 
given that the heliospheric/geocoronal solar-wind charge-exchange (SWCX) emission shows 
a hint of an Fe L-shell complex around ~0.8keV (Snowden et al. 2004; Fujimoto et al. 
2007; Smith et al. 2007). Therefore, the inferred intensities of O Ly(/3 + ^ + 6 + 6) should be 
considered to be upper limits. This model gives us statistically acceptable fits for all the cases 
examined. Figure M shows the best-fit model and the residuals. In the figure, left (labeled 
A) and right (labeled B) columns present results for the vpshock model with the riet-value 
of 2.8x10^^ cm~^sec~^ (see. Table [2]) and that with half the value, respectively. The best-fit 
EMs (not shown in Table H]) are obtained to be ~3xl0^^cm"^ for NE4, ~2.3xl0^'^ cm~^ for 
P27, ~1.4x 10^^ cm~^ for Rim2, respectively. 

In order to clarify improvements from the initial fit with the pure thermal emission 
model (Fig. [3] (a)), we show the NE4 spectrum with the initial (pure thermal) model in 
Fig. [7] (a) and the new (thermal plus possible CX lines) model in Fig. [7] (b), focusing on 
a small spectral region around the possible CX lines of interest (i.e., 0.68-0.76 keV). It is 
clear from this figure that the new model fits the data much better than the initial one. For 
quantitative evaluation of the goodness-of-fits, we perform Kolmogorov-Smirnov (K-S) test 
for the two fit results shown in Fig. [7] (a) and (b). We find that the probabilities that the 
data are drawn from the initial and new best-fit models are ~15% and ~80%, respectively. 
Therefore, the new model is much better than the initial one, although both models cannot 
be fully ruled out according to this K-S test. We also re- fit the data in a wide energy band 
using a maximum likelihood method, and compute the Bayesian Information Criterion (BIC) 
defined as + khi{n), where k is the number of free parameters and n is the number of 
spectral bins. The BIC-value is lower for the new model than the initial one, which also 
shows the validity for introducing possible CX lines. 

Table m summarizes the best-fit Gaussian parameters, from which the intensity ratios of 
O He(7 + 5 + e) / O Hea are calculated to be ~0.097-0.110 in NE4, ~0.057-0.061 in P27, and 
~0. 066-0. 090 in Rim2. These are much higher than that expected for the thermal emission 
from the Cygnus Loop (~0.02 is expected at kT^ ~ 0.2 keV), but are also larger than those 
expected in CX emission (the ratio is likely to be ~0.04; Beiersdorfer et al. 2003; Wargelin et 
al. 2005). These inconsistencies imply that our modeling is not perfect and requires further 
sophistication which is beyond the scope of this paper. Here, it may be worth noting that the 



- 14 - 



ratios derived in our primitive modeling are closer to CX emission than to thermal emission, 
lending support to the idea that CX emission is present. 

Using the line intensities listed in Table IH we next estimate ion number ratios in the 
plasma, as performed for SWCX emission (e.g., Kharchenko et al. 2003; Krasnopolsky 2004; 
Snowden et al. 2004). We limit our estimation to only O VIII, O IX, and Ne X, whose line 
intensities (K-shell lines from O VII, O VIII, and Ne IX after CX reactions) are relatively well 
constrained. From Table HJ the intensity ratios of O VIII lines / O VII lines and Ne IX lines 
/ O VII lines are calculated to be ~0.095-0.184 and ~0.079-0.086 in NE4, ~0.006-0.026 and 
~0.059-0.069 in P27, and ~0.052-0.156 and ~0.040-0.053 in Rim2, respectively. Above, 
we summed K-shell fluxes from all transitions in the same species, since the total number 
of K-shell photons reflects the number of CX reactions. To convert the flux ratios to ion 
number ratios, we use SWCX cross sections for hydrogen listed in Wargelin et al. (2004). 
The cross section ratios of O IX / O VIII and Ne X / O VIII are computed to be 1.66 and 
1.53, respectively. Dividing the flux ratios by the cross section ratios, the ion number ratios 
for O IX / O VIII and Ne X / O VIII are derived to be ~0.057-0.111 and ~0.052-0.056 in 
NE4, ~0.004-0.016 and ~0.039-0.045 in P27, and ~0.031-0.094 and ~0.026-0.035 in Rim2, 
respectively. 

We compare these ion number ratios with those expected in the thermal plasma of the 
Cygnus Loop's rim, i.e., fcre=0.2keV. We investigate three ionization timescales: riet = (0.5, 
1, 2) xlO^^ cm~^ sec. As described in the appendix, the ion number ratios for O IX / O VIII 
and Ne X / O VIII are derived to be <0.33 and 0.01-0.02, respectively (here, we adopt the 
Ne/0 abundance ratio of 2 times the solar value which is typically reported in the Cygnus 
Loop, e.g., Tsunemi et al. 2007). Although the error range is large for the former ratio, it 
is consistent with our estimates ranging from ~0.004 to ~0.111 using SWCX cross sections. 
On the other hand, the latter ratio is signiflcantly lower than our measurement by a factor 
of ~1.3-5.6. This discrepancy suggests either overestimates of CX lines from He-like Ne, 
underestimates of CX lines from H-like O, or both. The overall intensity of O VIII lines is 
also somewhat uncertain because of the difficulty in separating O Lya at 0.654 keV and O 
He/9 at 0.666 keV as well as the possible contribution of Fe L-shell lines around 0.8 keV. Based 
on the best-flt parameters in Tables [3] (total emission) andH] (CX emission), the fractions of 
CX emission in total emission are calculated to be ~40% for O Hea and ~50% for Ne Hea. 
A higher fraction of CX emission in Ne Hea than in O Hea is difficult to understand, since 
the fraction of H-like ions (which yield Hea CX emission) is higher for O than for Ne; at 
the electron temperature of ~0.2keV in thermal equilibrium, a little less than half of O is 
He-like, while less than 10% of Ne is H-like and the rest is He-like. Taking account of NEI 
effects would make things worse. Therefore, Ne Hea CX lines are likely to be overestimated 
by our primitive CX modeling. This problem could be caused by inaccurate modeling of the 
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thermal emission, because the hne ratios strongly depend on the plasma condition which 
is difficult to infer in the CX- contaminated spectra. In addition to the uncertainty in the 
spectral modeling, the CX cross sections might be different from those we employed. Given 
these complexities, the derivation of accurate ion number ratios from CX emission is left for 
future work. At the very least, we can argue here that the ion fractions inferred from both 
CX emission and the thermal emission model are consistent within an order of magnitude. 

4-3.2. CX Emission from SNRs 

Although it is challenging to distinguish between CX emission and thermal emission in 
our data, we have constructed a spectral model for the CX emission, assuming that spectral 
shape of the thermal emission is similar with azimuthal angle. We here discuss the inferred 
CX-to-thermal flux ratios (e.g., O Hea ~40%) in terms of theoretical expectations. 

Wise & Sarazin (1989) were the ffist to perform detailed calculations of CX-induced 
X-ray emission in a SNR. They found that CX emission generally contributes only 10"'^ 
to 10~^ of the collisionally excited lines in the entire SNR, even when the ambient gas is 
completely neutral. For comparison, we estimate the total CX-to-thermal flux ratio in the 
Cygnus Loop. We assume that the total CX flux is 10 times larger than that in Table H] 
for NE4, since the number of FOV where we see CX emission is about ten (see. Fig. El (a)). 
The total thermal flux from the remnant is likewise inferred from the integrated thermal 
flux in all the spectral extraction regions in Fig. [TJ The CX-to-thermal ratio in the entire 
SNR is then inferred to be on the order of 10~^, which is significantly higher than that in 
the model of Wise & Sarazin (1989). This discrepancy requires either a reduction of the 
CX contribution in our spectral modeling in Section 4.3.1, improvements of the theoretical 
calculation, or perhaps both. In fact, physical parameters such as the shock velocity used in 
the calculation by Wise & Sarazin (1989) are different from those in the Cygnus Loop. 

More recently, Lallement (2004b) investigated spatial variations of CX emission in a 
SNR. According to Fig. 3 in Lallement (2004b), CX emission can comprise a few 10% of 
the total X-ray (0.1-0.5 keV) emission within the radius of 0.98i?shock < R < -Rshock, which 
corresponds to the outermost rim region in our analysis. This theoretical prediction agrees 
reasonably well with our measurement of O Hea ~40%. Although Lallement's calculation 
is specifically intended for the LMC SNR DEM L71, we believe that it is appropriate for 
a qualitative comparison with the Cygnus Loop, since the basic parameters for both the 
DEM L71 SNR and the Cygnus Loop are similar. 

Lallement (2004b) also presented the expected EM of thermal emission at the equivalent 
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radius, i.e., the radius where the flux of CX emission becomes equal to that of the thermal 
emission. The expected EM is expressed as 2x10^'' ecos^^^ a x(T,a)~^ Tie Vioo cm~^, where 
e is the ratio between the CX probability and the coUisional ionization probability, 6 is an 
incidence angle (see. Fig. 2 in Lallement 2004b), a accounts for the hot ion content, i.e., the 
temperature and metallicity, x(T,a) is the thermal emissivity ratio for a given temperature T 
and the metal abundance a relative to a temperature 1x10^ K and the solar abundance, Uc is 
the ambient density, and Vioo is the relative velocity between neutrals and the hot gas in units 
of lOOkmsec"^. Following the discussion in Lallement (2004b), we employ e = 1, a = 1, and 
x(T,a)= 1. The equivalent radius of ~ 0.98i?shock, which is the inner edge of the outermost 
region in Fig. [1] gives us cos~^^ of ~5. The ambient density is taken to be unity from 
measured values ranging from 0.3 to 4cm~^ near NE4 (Raymond et al. 2003; Leahy 2003). 
The relative velocity can be considered as the shock velocity of ~300kmsec~^ (Ghavamian 
et al. 2001; Salvesen et al. 2009). Substituting these quantities in the equation above, we 
derive the expected EM of ~ 5 x 10^® cm~^. Since the EM observed in the outermost rim in 
NE4 is estimated to be ~1.4-3 xlO^^cm"^ (see. Section 4.2.1), the expected and observed 
EMs are consistent within a factor of ~2-6. Given the considerable uncertainties in the 
spectral modeling as well as the ambient density, we believe that the EMs are in reasonable 
agreement. 

We have found general agreements between our observational results and the theoretical 
predictions by Lallement (2004b). Nonetheless, further detailed calculation of the CX emis- 
sion in a SNR is encouraged, since Lallement 's calculation is somewhat simplified in that it 
does not consider the ionization structure of the post-shock region, and thus would lead to 
an overestimate of the CX contribution. Likewise, more sophisticated spectral modeling of 
CX emission is also desirable for future work. 



4.4. Need for High-Resolution Spectra with the Astro-H SXS 

While we have proposed CX emission for the origin of the "enhanced" abundance in 
the rim regions of the Cygnus Loop and have presented supporting observational evidence, 
solid evidence for the presence of CX emission is still missing. Therefore, at this point, we 
cannot fully exclude alternative explanations. For example, there could really be a bimodal 
distribution of abundances in the rim of the Cygnus Loop, in spite of the difficulties described 
in Section 1. 

The Astro-H SXS will hopefully provide definitive evidence for CX emission by resolv- 
ing the O Hea triplets, i.e., resonance, intercombination, and forbidden lines. In thermal 
emission, the ratio of forbidden to resonance line strengths is less than unity. In contrast. 
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for CX emission the ratio is expected to be about 10 (Kharchenko et al. 2003). Laboratory 
measurements of CX reactions do show stronger forbidden hues than resonance hues (Beiers- 
dorfer et ah 2003), ahhough the forbidden-to- resonance ratio is not as large as 10 (which 
might be related to the fact that the authors used CO2 rather than H as electron-donor 
species). Thus, we expect that forbidden lines will be enhanced in the CX-contaminated 
outermost rim regions compared with those in the rest of the rim. So far, the northern 
bright portion of the Cygnus Loop is the only location where the triplets have been resolved, 
by high spectral resolution Einstein Focal Plane Crystal Spectrometer (FPCS) observations 
(Vedder et al. 1986). The resonance line there is much stronger than the forbidden line, 
suggesting that CX emission is not significant. Since the FPCS FOV (3'x30') is located well 
outside (or inside with respect to the center of the Loop) our Suzaku FOV shown in Fig [H 
it is reasonable to consider that signatures of CX emission if present would be washed out 
by thermal emission. Spatial variations of the line ratio, which provide the key information 
about the presence of CX emission, will be measured with the Astro-H SXS. 

With the current CCD data, it might be possible to find a forbidden/resonance line-ratio 
variation by measuring the energy separation between O Hea and O Lya (e.g., Lallement 
2009). However, such energy separations are expected to be less than a few eV, and thus 
are too small to detect with nondispersive CCD spectra. As we show in Tables [3] and HI the 
fit line center energies in the Suzaku spectra are randomly scattered within the energy-scale 
calibration uncertainties (±5eV: Koyama et al. 2007). Moreover, O Lya is significantly 
contaminated by O He/3 in part of the Cygnus Loop, and its contribution is highly variable 
from location to location. This makes even more difficult any firm measurement of the 
separation between O Hea and O Lya. Also, dispersed spectra from slitless X-ray grating 
spectrometers are almost useless for the Cygnus Loop because of its large extent, although 
the XMM-Newton RGS might provide better spectra for some bright knotty features. Since 
the Astro-H SXS is a nondispersive microcalorimeter, it will be the first instrument to provide 
us with high-resolution spectra of sufficient quality to resolve individual lines of the triplets. 
The high-resolution spectra will also make possible identification of the individual lines 
around 0.7keV. Therefore, the Astro-H SXS is being eagerly awaited for further studies of 
the Cygnus Loop. 

In addition, the Astro-H SXS will help us derive much more accurate absolute metal 
abundances than those from nondispersive CCD spectra. This is because we can directly 
measure equivalent widths of lines below 1 keV with SXS spectra, which is not possible with 
CCD spectra. 

While measuring the absolute abundance with CCD spectra is quite difficult, it should 
be noted that relative metal abundances are already well constrained unless the spectra 
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are heavily contaminated by CX emission. Our analysis suggests that CX emission is not 
significant in the "normal" abundance regions. The relative abundances obtained there are 
close to the solar ratio except for O: the abundance of O is depleted by an additional factor 
of 2 relative to other elements (e.g., Miyata et al. 1994). As Miyata et al. (2008) pointed 
out, the depiction of the O abundance might be due to resonance-line scattering which is 
especially important for O K-shell lines. Therefore, in order to obtain accurate absolute and 
relative abundances in the rim of the Cygnus Loop, we will need to carefully consider effects 
of both CX emission and resonance-line scattering. 

5. Conclusion 

Using Suzaku and XMM-Newton data covering nearly the entire rim of the Cygnus Loop, 
we found that outermost rims in N-NE, SE, and NW regions show a strong line feature at 
~0.7keV. Its anomalously strong intensity and spatial localization in the outermost rim 
lead us to propose that it originates from CX processes, most likely between neutrals and 
H-like O ions. If so, this is the first detection of CX-induced X-ray emission associated with 
the forward shock of a SNR. The regions showing the ~0.7keV emission correspond to the 
"enhanced" abundance regions, and we suggest that the ~0.7keV emission and possibly 
other CX lines could enhance the apparent metal abundances there when we fit the CX- 
contaminated spectra with pure thermal emission models. The introduction of the CX 
emission potentially resolves the mysterious abundance "enhancements" in the outermost 
rim regions of the Cygnus Loop. On the other hand, the abundance depletion in overall rim 
regions still remains an open question. 
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6. Appendix 

To estimate tlie ion fractions, we first calculate ion fractions at collisional ionization 
equilibrium usiug sherpi Secoud, we calculate n.^T) (i.e.. a product of iou fractiou 
and emissivity) of some selected lines as a function of ngt, using the NEIline software 
provided by R. Smith (applications of this software can be found in literature, e.g., Sasaki 
et al. 2004; Katsuda & Tsunemi 2006). Above, nz{nct) /nz is the ionization fraction of 
the ionization species responsible for the line, and ei{T), a function of temperature, is the 
intrinsic emissivity of the line. Then, assuming that the ei{T) is constant with riet, we 
roughly estimate ion fractions at various riet-values. To evaluate systematic uncertainties 
associated with this method, we use not only Ka lines but also K/? lines. This method does 
not directly give ion fractions for fully ionized ions, since they do not radiate line emission 
by collisional excitation. We thus assume that fully ionized ions take up the balance, which 
is a good approximation in the plasma condition of interest. 



http://cxc.harvard.edu/ciao/sherpa 
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Table 1. Summary of observations. 
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Table 2. Spectral-fit parameters in the outermost regions of NE4 and NE2, using the 

absorbed vpshock model. 
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0.14±0.02 


EM"[xl0i^ cm-5].... 


n oc+0.20 

0.36I„„4 


2 q+0-5 


1.5±0.3 


Energy shifts [eV] 


+0.4 (FI), -1.7 (BI) 


+0.4 (FI), -1.7 (BI) 


-2.9 (FI), +0.2 (BI) 


x7d.o.f 


1219/551 


802/505 


466/410 



Note. — "EM denotes the emission measure (= J nefindi). Tlic values of abundances are multiples of 
the solar values (Anders & Grevesse 1989). The errors represent 90% confidence ranges. 
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Table 3. Gaussian parameters in the NE4, NE2, P27, P19, Rim2, and P24's outermost 
rims, using an absorbed phenomenological model consisting of an underlying 
bremsstrahlung plus 18 Gaussians. 



Line 


NE4'' 


NE2 


P27" 


P19 


Rim2'' 


P24 


C Lya 
Center (eV) . . . 
Width (eV) . . . 
Normalization'' 


357(357)13 
28±2 
7606+5S? 

— 405 


351(354)±5 
30±3 
3215+309 


349(351)±8 

32I5 
7998+?25„ 


349(345)±6 
32l^ 


349(356)±3 
30±2 
7074+5^6 


347(348)±6 
291-1 
5421+«6 

— 778 


C Ly/3 + N Hea 
Center (eV) . . . 
Width (eV) . . . 
N orm ciiliz ctt i on ^ 


438(438)±1 
1411 


437(441)±2 
14±3 
591+?? 


435(438)+4 
14l^ 

1779+}?? 
^' '="-116 


432(429)±2 

1411 
2239+ }?S 


435(442)±3 
111? 
958+?? 


433(434)+^ 

Q+6 
^-10 

822+25 


N Lya + N He/S 
Center (eV) . . . 
Width (cV) . . . 
Normalization^ 


499(499)12 
0= 
730±20 


497(501)1!^ 
0= 
177±12 


499(501)14 
0= 
481±38 


499(496)12 
0= 
686±39 


499(506)115 
0= 
281±23 


495(496)1-2 
0= 
257±25 


O Hca 
Center (eV) . . . 
Width (eV) . . . 
N orm ctliz ctt i on ^ 


565(565)±0 
12±0 
5547+5? 

— 55 


564(567)±1 
12±1 
1659+34 


564(567)±1 
13±1 
4763l?8g 


565(561)±1 

lll^ 
45591^^5 


564(571)±1 
12±1 
2804+?? 


563(564)±1 
I2II 
2243+?® 


O Lya 
Center (eV) . . . 
Width (eV) . . . 
Nor m 3,11 z ftt ion ^ 


654(653)= 
0= 
602±17 


651(655)= 
0= 
208+10 


653(655)= 
0= 
258±28 


653(649)= 
0= 
607±32 


653(660)= 
0= 
259+16 


651(652)= 
0= 
178±18 


O He/3 
Center (eV) . . . 
Width (cV) . . . 
Normalization^ 


666(665)^= 
367±16 


663(667)= 
0= 
82±9 


665(667)= 
0= 
444±26 


665(661)= 
0= 
437±29 


665(672)= 
0= 
130±14 


663(664)= 
0= 
194±16 


O He7+etc. 
Center (cV) . . . 
Width (cV) . . . 
Normalization'^ 


704(704)±2 
22±2 
341+11 

— 16 


708(711)+s 
48+i° 


703(706)±5 
24l^ 
239I27 


711(708)±6 
23±7 
2O6I29 


705(712)+* 
23l« 
145+;« 

— 15 


718(719)+9 
47+" 

—9 


Fe L (3d-i>2p) 
Center (cV) . . . 
Width (eV) . . . 
Normalization'' 


726(725)= 
0= 


723(727)= 
0= 


725(727)= 725(721)= 

0= 0= 
Linked to Fe L (3d->-2p) 


725(732)= 
0= 


723(724)= 
0= 


OLy/3 
Center (eV) . . . 
Width (cV) . . . 
Normalization'' 


775(774)= 
0= 
99±6 


772(776)= 
0= 
35±4 


774(776)= 
0= 
54±8 


774(770)= 
0= 
134itl2 


774(781)= 
0= 
38±5 


772(773)= 
0= 
42±6 


Fe L (3s->-2p) 
Center (eV) . . . 
Width (eV) . . . 
Normalization'' 


822 (822) ±1 
0= 
138±4 


819(823)12 
0= 
57±3 


819(821)1^ 
0= 
93±6 


822(818)13 
0= 
247±9 


822(829)1^ 
0= 
67±4 


819(820)±2 
0= 
64±4 


Ne Hea 
Center (eV) . . . 
Width (eV) . . . 
Normalization'' 


912(912)±1 
15±1 
3641^ 


909(912)11 
19±2 
891^ 


913(915)±1 
13±2 
282±8 


912(908)±1 
15±2 
337±10 


915(921)±1 
14±2 
140±4 


912(913)±2 
13li 
111±5 
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Table 3 — Continued 



Line 


NE4'' 


NE2 


P27" 


P19 


Rim2'' 


P24 


Ne hya 
Center (eV) . . . 
Width (cV) . . . 
Noiiiia.iiza.li(.)ii^' 


1022(1021)= 

1()±2 


1019(1023)= 

0= 

7±1 


1021(1023)= 

0= 
12±2 


1021(1017)= 

0= 
2()±:-! 


1021(1028)= 

0= 
7±1 


1019(1020)= 

0= 
10±2 


Ne He/3 
Center (cV) . . . 
Width (cV) . . . 
Nornicilizcitiori^ 


1074(1073)= 
0= 
19±1 


1071(1075)= 
0= 
5±1 


1073(1075)= 
0= 
15±2 


1073(1069)= 
0= 
16±3 


1073(1080)= 
0= 
7±1 


1071(1072)= 
0= 
6±1 


Ne He7 
Center (eV) . . . 
Width (eV) . . . 
Normalization'' 


1127(1126)= 
0= 
12±1 


1124(1128)= 
0= 
2±1 


1126(1128)= 
0= 
8±2 


1126(1122)= 
0= 
17±3 


1126(1133)= 
0= 
5+1 


1124(1125)= 
0= 
5±1 


Ne Ly/3 
Center (eV) . . . 
Width (eV) . . . 
Normalization'' 


1210(1209)= 
0= 
5±1 


1207(1211)= 
0= 
3±1 


1209(1211)= 
0= 
7±2 


1209(1205)= 
0= 


1209(1216)= 
0= 
4±1 


1207(1208)= 
0= 
4±1 


Mg Hea 
Center (eV) . . . 
Width (eV) . . . 
Normalization'' 


1341(1341)±4 
12 (<19) 
8.2±1.1 


1346(1350)±5 
(<31) 
2.9±1.0 


1354(1356)±11 
36li 

6 2+^* 
"•^-l.B 


1343(1339)±5 
6 (<22) 
9.9±1.8 


1350(1357)±7 

97+10 
^'-13 

4.811;? 


1346(1347)±8 
39l?« 
6.6±1.3 


Mg He/9 
Center (eV) . . . 
Width (eV) . . . 
Normalization'' 


1571(1571)±17 
0= 
1.2±0.8 


1511(1514)±38 
0= 

0.5 (<1.2) 


1522(1525)±381 
0= 
1.52ll;° 


1578(1574)±395 
0= 

1.2 (<2.4) 


1574(1581)±394 
0= 

0.8 (<1.6) 


1576(1577)±394 
0= 

0.9 (<1.8) 


Si Hea 
Center (eV) . . . 
Width (eV) . . . 
Normalization'' 


1833(1833)111 
0= 
1.3±0.9 


1886(1889)±22 

0= 

0.8 (<1.7) 


1839(1841)±57 

0= 

1.4 (<2.4) 


1865(1861)±18 
0= 
2.8l?;| 


1872(1879)±37 
0= 
1 5+1° 


1806(1808)±43 
0= 
l.litO.9 



Note. — ""Enhanced" -abundance regions. ''In units of photons cm~^ s~i sr~i. ''Fixed values. Line center energies include effects 
of energy shifts determined by the fittings. Values in parentheses arc derived from the Bl CCD. The errors (90% confidence level) arc 
estimated with the other Gaussian parameters and the bremsstrahlung parameters fixed at the best-fit values. These errors do not 
include systematic uncertainties, such as ±5eV for the line center energies (Koyama et al. 2007). 
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Table 4. Gaussian parameters for the possible CX emission in the NE4, P27, and Rim2's 
outermost rims, using an absorbed combination model consisting of the vpshock plus 22 

Gaussians. 



Line NE4 P27 Rim2 



C Lya 
Center (cV) . . . 
Width (eV) . . . 
Normalization" 


A 

343{343)±13 


B 


A 

345(347)±17 
2728tl«I| 


B 


A 

330(337)+" 

30lf 
6346l««9 


B 

313(320)+9^ 
29lf 
748511^16 


C Ly/3 + N Hca 
Center (cV) . . . 
Width (cV) . . . 
Normalization" 


436(436)^4 
20+1 
869lfo 


433(433)±3 
24+4 
922ll°4 


435(437)11 
20t* 
1082^202 


426(428)17 
22I9 
^^•^-192 


428(435)+6 
18+« 
85ll^f 


419(426)+6 

21II 
8l8t{ll 


N Lya + N He,9 
Center (cV) . . . 
Width (eV) . . . 
Normalization" 


500(500)±3 
0" 
160±21 


499(499) ±2 
0" 
2561^1 


510(512)+-"^ 
0" 
113±37 


503(506)+j2 
0'' 
11411 


505(512)±7 
0" 
llOlai 


507(514)+^ 
0'' 
1251^5 


O Hea 
Center (eV) . . . 
Width (eV) . . . 
Normalization" 


563(563)±1 


20581^1 


563(563)±1 


2307i«^ 


562(565)±1 

2262+1^? 


562(564)±1 


1926±125 


560(567)±1 


1281±71 


559(566)±2 


1054±75 


O Lya 
Center (eV) . . . 
Width (eV) . . . 
Normalization" 


654(653)'' 


654(653)'' 
0" 
3691S 


653(655)'' 

0" 
(<76) 


653(655)'' 
0" 

10 (<138) 


653(660)'' 
0" 

45 (<121) 


653(660)'' 
0" 
134lli 


O He/3 
Center (cV) . . . 
Width (eV) . . . 
Normalization" 


666(665)'' 
Of" 
2261^° 


666(665)'' 
0'' 
246l?| 


665(667)'' 
0" 
324l« 


665(667)'' 
0'' 

396l?!6 


665(672)'' 
0" 
I25I8I 


665(672)'' 
0'' 

86 (<133) 


O He7 
Center (eV) . . . 

Width (cV) . . . 
Nomiahzation" 


698(697)'' 

Of- 

70 (<102) 


698(697)'' 

0" 
TOtll 


697(699)'' 

0" 

39 (<83) 


697(699)'' 

0" 
(<63) 


697(704)'' 

0" 
(<41) 


697(704)'' 

0" 
(<42) 


O He<5 
Center (eV) . . . 
Width (eV) . . . 
Normalization" 


713(712)'' 

(<iio) 


713(712)'' 

0'' 
(<67) 


712(714)'' 
0'' 

(<109) 


712(714)'' 

0'' 
(<87) 


712(719)'' 

0'' 
(<59) 


712(719)'' 

0'' 
(<58) 


O Hee 
Center (eV) . . . 
Width (eV) . . . 
Normalization" 


723(722)'' 
0'' 
1301^62 


723(722)'' 
O*" 
183lfo 


722(724)'' 

0" 
91+31 


722(724)'' 
0'' 
118+11 


722(729)'' 
0" 

°°-50 


722(729)'' 
0'' 
951^1 


O Ly/3 
Center (eV) . . . 
Width (cV) . . . 
NormaUzation" 


775(774)'' 

Of- 

40l'^ 


775(774)'' 

0'' 
76±10 


774(776)'' 

0'' 


774(776)'' 
0'' 
29±14 


774(781)'' 

Of- 

18±8 


774(781)'' 
0'' 
2911° 


O Ly-/ 
Center (eV) . . . 
Width (eV) . . . 
Normalization" 


817(816)'' 
0'' 
(<5) 


817(816)'' 

0'' 
of5+16 
•30_24 


816(818)'' 

0'' 
(<10) 


816(818)'' 
0'' 

11 (<32) 


816(823)'' 

0'' 
(<11) 


816(823)'' 

0'' 
9 (<23) 



Table 4 — Continued 



Line NE4 P27 Rim2 



O LyS 



Center (eV) . . . 
Width (eV) . . . 


837(836)'' 
0" 


837(836)'' 
0" 


836(838)'' 
0" 


836(838)'' 
0" 


836(843)'' 
0" 


836(843)'' 
0" 


Normalization" 


«11) 


14 (<52) 


(<12) 


14 (<28) 


6 (<18) 


11 (<66) 


O Lye 














Center (eV) . . . 
Width (eV) . . . 


849(848)'' 
0" 


849(848)'' 
0" 


848(850)'' 
0" 


848(850)'' 
0" 


848(855)'' 
0" 


848(855)'' 
0" 


Normalization"^ 


16lti 


20 (<35) 


5 (<12) 


(<17) 


8 (<17) 


10 (<23) 


Ne Hea 














Center (eV) . . . 
Width (eV) . . . 


913(913)±1 
13±2 


914(914)±1 
14±2 


914(916)±2 
6l^ 


914(916)+! 


918(925)+2 


918(925)+2 
10+5 


Normahzation" 


181±5 


221±5 


150±8 


156t« 


57±4 


61±5 


Ne hya 














Center (eV) . . . 
Width (eV) . . . 


1022(1021)'' 
0'' 


1022(1021)'' 
0'' 


1021(1023)'' 
0'' 


1021(1023)'' 
0'' 


1021(1028)'' 
0'' 


1021(1028)'' 
O*" 


Normalization" 


3±2 


10±2 


2 (<5) 


4±3 


(<2) 


2±1 


Ne He/3 














Center (eV) . . . 
Width (eV) . . . 


1074(1073)'' 
0" 


1074(1073)'' 
0" 


1073(1075)'' 
0" 


1073(1075)'' 
0" 


1073(1080)'' 
0" 


1073(1080)'' 
0" 


Normalization" 


10±2 


13±2 


7+2 
'-3 


8±3 


2t? 


3±2 


Ne He7 














Center (eV) . . . 
Width (eV) . . . 


1127(1126)'' 
0'' 


1127(1126)'' 
0'' 


1126(1128)'' 
0" 


1126(1128)'' 
0" 


1126(1133)'' 
0" 


1126(1133)'' 
0" 


Normalization" 


3±? 


6l? 


4±2 


5±2 


(<4) 


1 «3) 


Ne He<5 














Center (eV) . . . 
Width (eV) . . . 


1150(1149)'' 
0" 


1150(1149)'' 
0" 


1149(1151)'' 
0" 


1149(1151)'' 
0" 


1149(1156)'' 
0" 


1149(1156)'' 
0" 


Normalization" 


2±1 


2±1 


(<2) 


(<1) 


(<2) 


1 (<4) 


Ne Ly/3 














Center (eV) . . . 
Width (eV) . . . 


1210(1209)'' 
0'' 


1210(1209)'' 
0'' 


1209(1211)'' 
0" 


1209(1211)'' 
0" 


1209(1216)'' 
0" 


1209(1216)'' 
0" 


Normalization" 


2±1 


4±1 


1 «2) 


2±1 


0«1) 


0«1) 


Mg HcQ 














Center (eV) . . . 


1343(1343)lg 


1342(1342)tg 






1350(1357)±10 


1346(1353)±15 


Width (eV) . . . 


(<16) 


13 (<24) 


(<103) 


(<75) 


(<1796) 


17 (<42) 


Normalization" 


3.1±1.1 


5.5±1.2 


1.3 (<2.9) 


2.0l« 


-■-■3-1.0 


^■-'--1.0 


Mg He/3 














Center (eV) . . . 
Width (eV) . . . 


1575(1575)±22 
0" 


1570(1570)+^^ 
0" 


1538(1540)+|g 
0" 


1537(1540)+|o 
0" 


0" 


0" 


Normalization" 


0.8 (<1.6) 


1.2±0.8 


0.7 (<1.7) 


0.8 (<1.8) 


0.1 (<0.9) 


0.2 (<1.0) 


Si Hea 














Center (eV) . . . 
Width (eV) . . . 


1835(1835)±15 
0" 


1834(1834)±15 
0" 


0" 


0" 


1866(1873)±36 
0" 


1864(1871)±34 
0" 


Normalization" 


U.M_(j 7 


1 0+^'" 


0.8 (<1.9) 


0.8 (<1.9) 


1.0±0.8 


l.litO.8 



Note. — Columns A and B are responsible for the thermal emission models with net = 2.8 X 10^^ cm sec ^ and net = 
1.4 X 10^^ cm~^sec~^, respectively. Other notes can be found in TablefS] 
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Fig. 1. — Logarithmically-scaled ROSAT }iIU. image of the entire Cygnus Loop. The image 
is binned by 5" and has been smoothed by a Gaussian kernel of a = 15". The FOV analyzed 
in this paper are shown as boxes for Suzaku and a circle for XMM-Newton. For convenience, 
we label Suzaku FOV according to the object name in the data file. Spectral extraction 
regions are indicated as narrow lines along the shock front. XIS spectra extracted from 
the outermost regions in the FOV in cyan ("enhanced" abundances) and green ("normal" 
abundances) are shown in Fig. |6l 
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Ovii 




Energy (keV) 

Fig. 2. — Two FI (XISO+2+3) spectra extracted from the outermost rim regions in the 
NE2 (black) and the NE4 (red). The NE2 spectrum is muhiphed by 2.33 so that the count 
rate in the continuum band (1.4-1.5 keV and 1.65-1.75 keV) is equahzed to that in the NE4 
spectrum. The O and Ne enhancement in the NE4 spectrum is readily seen. 




0.5 12 0.5 1 2 

Energy (keV) Energy (keV) 



Fig. 3. — (a) XIS spectrum extracted from the NE4's outermost rim along with the best-fit 
single temperature vpshock model. Black and red are the FI (XISO-l-2-1-3) and BI (XISl), 
respectively. Contributions of the X-ray backgrounds (see, text) are separately shown only 
for the BI data. The residuals are shown in the lower panel, (b): Same as left but for NE2. 
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0.5 1 2 

Energy (keV) 

Fig. 4. — (a) NXB-subtracted FI (XISO+2+3) spectrum from the NE4's outermost rim (FI 
only) along with the best-fit model consisting of a bremsstrahlung plus 18 Gaussians. Line 
identifications are shown on top of the panel. Green curves show individual contributions of 
the X-ray background, (b) Same as above but for the NE2's outermost rim. 
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5 a.2 0.4 Oe 0£ 1 1.2 14 1£ 1£ 2 

Fig. 5. — (a) Line-ratio map of the ~0.7keV emission to the Fe L-shell complex at ~0.82 keV. 
The X-ray boundary of the Cygnus Loop is shown as a white hne. Arrows indicate the 
regions where we show the XIS spectra in Fig. |6l (b) Logarithmically-scaled radio (325 
MHz) surface brightness map obtained by the WENSS team (Rengelink et al. 1997), with 
the X-ray boundary of the Cygnus Loop as in Fig. [5] (a). 
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Fig. 6. — Top panels: NE4's outermost rim spectra (black and red for the FI and BI data, 
respectively) along with the best-fit model consisting of an assumed vpshock component 
(blue), 22 Gaussians (black), and the X-ray background (green). Contributions of these 
components are separately shown only for the BI data. The vpshock component in the left 
panel has the same spectral parameters derived in the NE2's outermost rim (see. Table [2]), 
whereas that in the right panel has half the value for the ionization timescale. Middle and 
bottom panels: Same as top panels but for P27 (middle) and Rim2 (bottom). 
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Fig. 7. — (a) The NE4's outermost rim spectrum (only for FI) along with the best-fit model 
consisting of pure thermal emission (see, Table [2]). This is identical to that in Fig. [3] but 
focusing on a small spectral region of 0.68-0.76 keV. (b) Same as Fig. [7] (a) but the model 
includes the possible CX lines. 
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